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ral  angle  (less  than  60  degrees)  there  is  an  equilibrium  porosity  in  a,poLy- 
crystalline  material;  for  sucli  a  material  there  is  no  barrier  to  nucleation  of 
cavities  in  dense  material.  At  intermediate  dihedral  angles,  pores  with  few 
sides  would  shrink  and  disappear,  whereas  those  with  many  sides  would  initiall 
shrink  but  would  reach  an  equilibrium  size  at  which  dens i f ication  would  stop. 
In  this  case  grain  growth  of  the  large  pores  would  reduce  the  number  of  sides 
for  the  pore  and  promote  further  densification.  For  such  materials  there 
would  be  a  small  barrier  to  nucleation.^ — 
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1.  Introduction 


This  report  describes  progress  which  has  been  made  during  he  xast 
year  on  oxidation  and  sintering  of  non-oxide  ceramics.  A  number  of  personnel 
activities  occurred  in  this  period.  Ms.  Monaghan  who  joined  us  in  August, 

1979  from  Birmingham  University  to  work  on  the  oxidation  of  Si^Nj  departed 
in  the  Summer  of  1980.  Ms.  Mieskowski  joined  us  in  January  1980  from  Penn 
State  University  via  Union  Carbide  and  continues  to  work  on  oxidation  of 
SiC  and  will  soon  write  her  M.S.  thesis.  Prof.  Cannon  spent  January  to  June 

1980  on  sabbatical  leave  from  MIT  to  do  theoretical  work  on  sintering  of 
non-oxide  ceramics  in  cooperation  not  only  with  us,  but  also  with  Prof. 
Moller  who  had  been  on  leave  from  the  University  of  Gottingen  working  on 
grain  boundaries  and  sintering  of  silicon.  Prof.  Moller  has  since  left  but 
stays  in  contact.  Dr.  Tighe  has  continued  to  visit  from  NBS  to  perform 
oxidation  experiments  on  Si^  using  the  environmental  cell  of  the  high 
voltage  electron  microscope  (HVEf)  . 

2.  Oxidation  of  SiC 

During  the  past  20  years,  a  number  of  investigations  concerning  the 
mechanism  of  SiC  oxidation  have  been  conducted.  In  general,  these  studies 
have  supported  the  idea  that  oxidation  occurs  by  the  inward  diffusion  of 
some  species  of  oxygen  through  the  SiO^scale.  The  initial  work  on  this 
contract  raised  the  question  of  whether  this  concept  was  in  error  based  on 
observations  of  spherulitic  crystallization  of  the  amorphous  SiO^scale 
(See  our  publication  in  Appendix  1) . 

Spherulites  in  oxide  scales  <1  pm  in  thickness  are  actually  disc-shaped 
rather  than  spheriodal.  They  have  been  observed  on  all  of  the  samples 
tested  to  date  which  have  included  CVD  and  Acheson-type  single  crystals  and 
a  and  B  polycrystalline  SiC.  However,  growth  of  additional  layers  of 
spherulites  on  top  of  previously  existing  ones  was  not  apparent  in  any  of 
these  experiments. 
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A  tracer  study  involving  CVD  single  crystals  oxidized  both  in  ^0  and 

18 

0  was  performed  but  was  unsuccessful  in  indicating  the  direction  of  the 
diffusing  species.  The  samples  had  been  oxidized  as  received  instead  of  being 
polished  sufficiently  to  assure  a  flat  surface:  data  interpretation  was 
impossible  because  of  the  crystal  growth  steps.  The  experiment  will  be 
repeated  but  with  specimens  of  a  and  B  polycrystalline  SiC  which  have  already 
been  polished  and  oxidized  in  ^0. 

An  experiment  is  also  being  planned  which  will  provide  information  on 
the  effect  of  cleanliness  of  spherulitic  growth.  All  specimens  to  date 
have  been  cleaned  with  dilute  HF,  distilled  1^0,  and  ethanol  under  open 
laboratory  conditions  and  oxidized  in  furnaces  which  are  used  for  various 
types  of  controlled  atmospheric  work.  Selected  specimens  will  be  prepared 
in  the  "clean  room"  used  for  microelectronics  research  and  then  oxidized  in 


an  adjacent  area.  These  will  then  be  examined  and  compared  to  previously 
oxidized  samples. 

The  following  table  is  a  summary  of  recent  oxidations  and  microscopy: 


Single  crystals 


Temp.  Time 


1300°C  1,2,5,10,20  ambient  CVD,  Acheson  optical 

40,60  hrs.  (also  single 

crystal  Si) 


1350°C  1,2,3,5,10,  1  atm  CVD,  Acheson 

20  hrs. 


1400°C  88  hrs. 


ambient  CVD,  Acheson 


,«*  orr.es  or 

rones  or  waiintt  « at  ■* 

omriwti-  u  j 


i 


Poly crystal ilne  SiC 


Temp. 

Time 

Thickness 

°2 

SiC  type 

microscopy 

1200°C 

* 

0-0.25,0.50 

1.0pm 

1  atm 

a,B 

optical. 

SB1 

1300 °C 

* 

o<0.25,0.50 

1.0pm 

II 

II 

optical. 

SEM,  TEM 

1350°C 

28,32, 

128  hrs . 

* 

ambient 

1? 

optical 

1400°C 

* 

0-0.25,0.50 

1.0 

1  atm 

II 

optical , 

SB-1 

*  not  a 

primary  consideration 

The  SEM  examination  of  oxidized  a  and  3  polycrystalline  SiC  revealed 
bubbles  which  may  be  attributed  to  CO,  and  perhaps  SiO,  formation  at  the 
substrate/scale  interface.  EDAX  will  be  used  to  check  bubble  composition. 

TEM  studies  have  centered  on  spherulite  morphology  and  the  associated 
defect  structure.  The  spherlite  are  found  to  be  heavily  faulted  and  twinned, 
due  to  the  B  -*■  a  cristobalite  transformation.  EDAX  will  be  used  to  determine 
the  extent  of  impurity  atom/ion  migration  to  the  spherulite  grain  boundaries. 

3.  Oxidation  of  Si^N^ 

The  basis  of  this  work  is  given  in  Appendix  2,  which  is  a  reprint  of  our 
work  published  in  Electron  Microscopy  1980.  We  are  continuing  to  co-operate 
with  Dr.  Tighe  to  use  the  HVB1  environmental  cell  and  analytical  electron 
microscopy  techniques  to  understand  the  catastrophic  oxidation  of  ¥£0^- doped 
SijN^.  This  work  has  been  quiescent  since  Ms.  Monaghan  left  but  will  pick  up 
now  that  H.  C.  Liu  (who  has  just  finished  his  Ph.D.  with  Prof.  Mitchell)  has 
agreed  to  stay  on  and  dc  this  research. 

4.  Sintering  of  Covalent  Solids. 

Two  aspects  of  sintering  have  been  considered  from  a  theoretical  point 
of  view.  One  is  on  the  structure  of  grain  boundaries  in  covalent 
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solids  using  the  coincidence  site  lattice  model  as  applied  to  Si.  This 
work  was  done  by  Dr.  Moller  during  his  stay  here  and  the  published  paner 
is  given  in  Appendix  3.  The  other  area  is  on  the  influence  of  pore  curvature, 
dihedral  angle,  applied  stress  and  internal  pore  pressure  on  densification, 
in  an  attempt  to  understand  the  poor  sinterability  of  covalent  solids.  This 
work  was  undertaken  in  co-operation  with  Professor  Cannon  during  his 
Sabbatical  leave  here  last  year  and  the  first  draft  of  a  paper  is  given  in 
Appendix  4. 

5.  Other  Research. 

Several  other  papers  have  been  published  over  the  past  year  resulting 
from  AFOSR- sponsored  research.  These  include: 

(a)  A.H.Heuer,  L.U.Ogbuji  and  T.E. Mitchell,  "High  Resolution  Studies  of 
the  Final  Stages  of  the  8  ■*  a  Transformation  in  Polycrystalline  SiC", 
in  "Electron  Microscopy  and  Analysis  1979",  (Inst,  of  Physics,  London), 
p.  453  (1979). 

(b)  L.U.Ogbuji,  T.E. Mitchell  and  A.H.Heuer,  "The  6  -*•  a  Transformation  in 
Polycrystalline  SiC:  III,  The  Thickening  of  a  Plates",  J.  Amer. 

Ceram.  Soc.,  6£,  91  (1981). 

(c)  L.U.Ogbuji,  T.E. Mitchell,  A.H.Heuer  and  S.  Shinozaki,  "The  8  -*•  a 
Transformation  in  Polycrystalline  SiC:  IV,  A  Comparison  of  Conventionally 
Sintered,  Hot-Pressed,  Reaction-Sintered,  and  Chemically  Vapor-Deposited 
Samples",  J.  Amer.  Ceram.  Soc.,  6£,  100  (1981). 


- A - 


354 


Reprinted  from  the  Journal  of  The  American  Ceramic  Society,  Vol.  63,  No  5-6  Mey-June,  i960 
Copyright  1980  by  The  American  Ceramic  Society 


Appendix  1 


Journal  of  The  American  Ceramic  Society —Discussions  and  Notes  Vol.  63,  No.  5-6 


The  Mkrostructure  of  Oxide  Scales  on  Oxidized  Si  and 
SiC  Single  Crystals 

A  H  HEUER.*  L.  U.  OGBUJ1.*  and  T.  E.  MITCHELL* 


The  first  oxidation  product  of  both  Si  and  SiC  is  amorphous  SiOj. 
which  devitrifies  to  cristobalite  during  growth  of  the  scale.  The 
present  note  describes  preliminary  studies  of  oxidation  of  single 
crystals  in  unmodified  laboratory  ambients,  at  1 100°  to  1250*  for  Si 
and  1350°  to  1500°C  for  SiC;  the  cristobalite  invariably  crystallizes 
with  a  spherulitic  morphology.  The  initially  formed  sphemlites 
provide  a  convenient  marker  for  ascertaining  whether  oxidation 
occurs  at  the  substrate-SiO,  or  the  SiOj-gas  interface.  The  Si  and 
SiC  single  crystals  differ  most  significantly  in  this  regard:  in  SiC, 
after  the  first  layer  of  spherulites  forms ,  subsequent  spherulites  form 
on  top  of  the  previous  layer;  in  Si  the  initial  appearance  of  the  first 
layer  remains  unchanged  with  further  oxidation,  suggesting  that 
subsequent  oxide  forms  beneath  the  first  layer. 

Various  grades  of  a-SiC  single  crystals  were  investigated,  rang¬ 
ing  from  a  relatively  impure,  commercial-grade  material,*  through 
variously  doped  samples  to  colorless,  semiconductor-grade  materi¬ 
al.*  The  B-doped  Si  samples  were  also  of  semiconductor  grade  * 
Both  materials  were  sectioned  along  their  close-packed  planes 
((0001 )  for  SiC.  {111}  for  Si),  polished,  and  oxidized  for  various 
times  in  the  laboratory  ambient.  The  results  shown  here  are  from 
SiC  samples  oxidized  at  1400°  and  Si  samples  oxidized  at  1250°C. 
but  are  typical  of  all  the  samples  investigated.  The  structure  of  the 
oxide  scale  was  identified  by  X-ray  and  electron  diffraction  and  its 
morphology  was  determined  by  optical  and  transmission  electron 
microscopy  (TEM).  All  scales  studied  gave  rise  to  interference 
colors,  from  which  we  infer  a  scale  thickness  £  1  gm;  the  spheru- 
lites  are  thus  actually  disks  which  are  quite  thin  compared  to  their 
diameters  ( =»  100  pim). 
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Figure  1  is  an  optical  micrograph  of  cristobalite  spherulites  on 
SiC  after  oxidation  for  1 2  h  at  1 400°C,  showing  that  smaller  spheru¬ 
lites  grow  on  top  of  the  larger  ones.  The  sequence  of  such  growth  is 
illustrated  in  Fig.  2.  Figure  2(4)  shows  the  initial  layer  of  spherulites 
formed  after  1  h  at  1400°C;  additional  oxidation  for  1 1  h  produced 
the  second  generation  of  finer  spherulites  evident  in  Fig.  2(B)-  The 
result  of  a  similar,  sequential  oxidation  experiment  on  Si  at  1 250*0 
is  shown  in  Fig.  3,  where  Fig.  304 )  shows  the  structure  after  20  h, 
whereas  Fig.  3(B)  shows  the  same  region  of  scale  after  52  h.  The 
gross  surface  appearance  of  the  oxide  scale  on  Si  did  not  change 
with  oxidation  time;  TEM  investigation  showed  that  the  change  in 
the  apparent  “texture”  of  the  spherulites  on  Si  involved  reduction 
of  the  amount  of  residual  glass  within  the  spherulites  and  coarsening 
of  their  constituent  cristobalite  microcrystals. 

Two  basic  kinetic  questions  with  respect  to  the  mechanism  of 
oxidation  of  these  materials  are  of  interest:  first,  whether  oxidation 
occurs  by  inward  diffusion  of  oxygen  and  reaction  at  the  substrate/ 
SiOj  interface  (with  subsequent  outward  diffusion  of  CO  or  some 
other  carbonaceous  species  in  the  case  of  SiC)  or  by  outward 
diffusion  of  Si  in  the  case  of  Si  and  Si  and  C  (or  SiO  and  CO)  in  the 
case  of  SiC  and  reaction  at  the  Si02/air  interface  and,  second, 
whether  the  rate-controlling  process  is  inward  diffusion,  outward 
diffusion,  or  an  interface-reaction  mechanism.  The  mtcrostructural 
evolution  of  the  oxide  scale  on  Si  suggests  that  oxygen  diffuses 
through  the  oxide  scale  to  react  at  the  Si/Si02  interface,  a  point  of 
view  consistent  with  the  previous  literature  on  oxidation  of  Si.1'1 


F!|.  I .  Optical  micrograph  of  apherulitk  cristobalite  scale  on  SiC  oxidized 
for  12  h  at  I4OTC 


Fig.  2.  Optical  micrographs  of  spherulitic  cristobalite  Kale  of  SiC  (A) 
oxidized  for  1  h  at  l40<rC  and  (B)  further  oxidized  for  1 1  h  at  1400* . 
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Our  results  on  SiC  suggest  strongly,  however,  that  oxidation  of  this 
material  occurs  at  the  Si02/air  interface,  not  at  the  Si02/SiC  inter¬ 
face.  This  outw  ard-grow  ing  scale  requires  that  outward  transport  of 
Si  and  C  in  some  form  (SiC  "molecules."  SiO.  CO,  etc.)  be  the 
rate-controlling  process  during  oxidation  of  SiC.  This  point  of  view 
is  contrary  to  many  assertions  in  the  literature  which  imply  that  the 
inward  diffusion  of  oxygen  is  the  rate -controlling  process  That  the 
actual  rate -controlling  process  for  oxidation  of  SiC  is  controversial 
can  be  seen  from  the  following 

to)  Motzfeldt3  examined  the  older  data  on  the  oxidation  of  Si 
and  SiC  and  concluded  that,  in  both  cases,  inward  diffusion  of  02  is 
rate -controlling 

(/> )  Jorgensen  et  al.  *  discounted  the  idea  of  molecular  diffusion 
of  02.  CO.  or  C02  and  suggested  that  the  diffusing  species  must  be 
O  or  St  ions. 

(<  I  Fitzer  and  Ebi5  found  that  the  Arrhenius  plot  of  parabolic 
oxidation  rate  had  three  stages  as  a  function  of  reciprocal  tempera¬ 
ture  and  proposed  that  the  rate-controlling  mechanism  changes  from 
Oa  diffusion  below  =  I400°C  to  CO  diffusion  at  higher  tempera¬ 
tures  . 

id  I  Singhal6  found  a  very  high  activation  energy  for  the 
parabolic  oxidation  rate  (480  kJ/mol  vs  =60  to  280  kJ'mol  found  by 
previous  investigators)  and  concluded  that  desorption  of  CO  from 
the  StC-St02  interface  must  be  rate-controlling 


None  of  the  previous  investigators  systematically  considered  the 
effects  of  impurities  on  the  oxidation  kinetics,  although  we  have 
shown7  that  high-temperature  annealing  of  sintered  SiC  improves 
the  oxidation  resistance  (presumably  by  evaporation  of  the  B  added 
as  a  sintering  aid)  and  it  is  common  knowledge  that  the  oxidation 
resistance  of  hot-pressed  SiC  depends  on  the  type  and  amount  o! 
densification  aid  added  Furthermore,  there  has  been  inadequate 
consideration  given  to  the  effect  of  devitrification  on  either  the 
kinetics  or  mechanism  of  oxidation.  It  is  thus  apparent  that  the 
oxidation  of  SiC  is  still  far  from  understood,  but  it  appears  from  the 
present  work  that  outw  ard  diffusion  of  Si  andC  is  important  It  is  not 
certain,  however,  whether  Si  or  C  (or  SiO  or  CO)  is  the  rate- 
controlling  species. 

Note  Added  is  Prooe:  Costello  ei  ol .  *  recently  reported  that 
platinum  marker  experiments  showed  that  the  oxide  scale  on  SiC 
grew  by  inward  diffusion  of  oxygen  and  that  the  spherulites  were 
nucleated  at  the  SiC  surface.  Their  work  thus  suggests  (hat  the 
second  generation  of  spherulites  observed  in  our  experiments  might 
have  occurred  by  the  crystallization  of  vitreous  silica  remaining  on 
top  of  the  first  generation  of  spherulites  It  is  clear  that  the  oxidation 
of  SiC  is  far  from  understood;  experiments  are  in  progress  using 
isotopic  tracers  to  attempt  to  resolve  some  of  these  points 
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Hot-pressed  Si  3N4+8-I3W2O  3  compacts,  vhich  exhibited  catastrophic  oxidation  in  air 
at  7^0 °C  but  only  passive  oxidation  at  1380°C,  were  studied  by  HVEM  using  an  environ¬ 
mental  cell.  The  in  situ  experiments  were  started  as  part  of  a  larger  program  to 
establish  a  model  for  this  catastrophic  oxidation  process.  The  low  temperatures 
required  for  initiation  of  oxidation  are  within  the  capability  of  the  hot  stage  and 
the  reaction  products  are  sufficiently  complex  to  require  analysis  by  electron  dif¬ 
fraction.  This  yttria-doped  silicon  nitride  ceramic  is  of  interest  for  use  in  cerami: 
turbine  engine  components  because  of  its  good  mechanical  properties  at  temperatures 
above  120C°C  [l].  However,  until  the  cause  of  the  low  temperature  catastrchpic  oxida¬ 
tion  is  identified  and  corrected  the  material  cannot  be  manujaetured  successfully  for 
these  components.  The  lew  temperature  oxidation  is  similar  in  appearance  to  the  so- 
called  "pest  condition"  and  occurs  in  some  billets  and  not  in  others  which  were 
manufactured  under  the  same  apparent  commercial,  powder  processing  and  hot  pressing 
conditions.  Lange  ei  al.  [2]  found  that  billets  they  prepared  in  the  Si  3U4 — S  i  2  *■'  2  0-V  2 "  1 
compatibility  triangle  of  the  phase  equilibrium  diagram,  shown  in  Fig.  1,  were  oxida¬ 
tion  resistant;  and,  that  billets  with  compositions  outside  this  triangle  were  unstable 
in  oxidizing  conditions. 

In  this  paper,  the  material  designated  -A-  exhibited  only  passive  oxidation  during 
heating  in  air  from  600°C  to  1L0Q°C;  while  the  material  designated  -B-  exhibited 
catastrophic  oxidation  during  heating  in  air  at  7**0°C.  Both  billets  are  from  the  same 
manufacturer  [3].  A  bar  of  material  -B-  when  heated  in  a  gradient  furnace  from  735°C 
to  1360°C  for  20  hours  oxidized  catastrophically  at  the  low  temperature  end,  and 
passively  at  the  high  temperature  end.  The  major  crystalline  phases  in  both  materials 
were  identified  by  powder  x-ray  diffraction  as  BS^Ni,,  Yx  qS!  7O2  3Nt,  (H  phase)  and  WSi2 
(there  were  some  unidentified  lines).  The  tungsten  phase  (3-k%)  is  a  result  of 
contamination  from  WC  balls  used  in  ball-milling  of  the  silicon  nitride  powder,  and 
thus  is  a  common  impurity  phase  in  hot-pressed  S^Nt,.  The  WSi2  phase  was  not  present 
in  the  catastrophically  oxidized  portion  of  the  material  B,  but  was  present  in  the 
portion  oxidized  at  1380°C.  Although  the  compositions  of  A  and  B  were  similar,  there 
was  considerable  phase  segregation  and  inhomogeneity  in  material  B  as  seen  in  the 
light  microscope.  3oth  materials  appear  to  have  compositions  outside  the  compati¬ 
bility  triangle  discussed  by  Lange  rt  al.  [2]. 

The  Swann  environmental  cell  in  the  CWRU  650kV  electron  microscope  was  used  with  a 
platinum  strip  heater,  in  flowing  oxygen.  Specimens  of  materials  -A-  and  -B-  were 
prepared  as  ion-thinned  discs  with  one  flat  side  and  one  dished  side.  This  configura¬ 
tion  gave  good  contact  with  the  grid  heater;  and,  resulted  in  better  heat  conduction 
to  the  specimen. 

The  microstructure  of  the  yttria  -A-  is  shown  in  Fig.  2,  and  is  similar  to  that 
found  in  other  billets  [1,5].  The  large  tabular  grains  are  Si3Nc,,  the  triangular 
shaped  phases  surrounding  the  grain  corners  are  the  ¥2^3  +  phases  and  the  dark 

particles  along  the  grain  boundaries  are  the  WSij  and/or  WC  phases.  During  hot- 
pressing  at  1750-l800°C,  the  yttrium  oxynitride  phases  can  melt  and  form  a  series  of 
solid  solutions  with  the  impurity  phases  and  with  the  silica  in  the  starting  powder 
(5).  Both  crystalline  and  non-crystalline  phases  of  variable  composition  can  fill  the 
interstices  around  Si3.N1,  grains. 
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Crecinens  were  heated  gradually  to  ?50°C  in  uQ  torr  of  r'2  .  Under  these  cor. ;  • 
oxide  platelets  were  seen  to  nucleate  on  the  surfaces  of  the  grains,  ar  . 

holes  and  at  the  yttria  phases  as  in  Fig.  3.  The  platelets,  which  were  random^, 
oriented  and  incoherent  with  the  matrix,  gave  spotty  ring  patterns  both  at  ter.:-  •  .* 
ar.i  after  cooling  to  room  temperature.  The  platelets  in  Fig.  were  identifier  -  - 
these  ring  patterns  as  a  mixture  of  Si02  and  C^I^O.  During  the  oxidation  cmci  . . 
able  drift  and  fle*ure  ^ f  the  specimens  occurred,  as  is  evident  in  the  figures. 

After  a  few  minutes  of  heating,  cracks  appeared  along  some  grain  boundaries;  : 

as  oxidation  proceeded  a  film  was  observed  to  grow  by  the  merging  of  nucleating 
clusters.  These  clusters  changed  diffraction  contrast  as  they  merged.  It  is  as. 
at  this  time,  that  this  process  represented  the  gradual  reaction  and  volaticat ior. 
the  tungsten  rich  phase.  The  film  thickened  and  the  crack  widened  to  form,  a 
hole  as  seen  in  Fig.  Iron  rich  and  tungsten  depleted  phases  l  .d  reacted  to... 

foil  edges  were  detected  with  an  EDX  system  on  a  200  KV  instrument  [7]. 

The  yttria-siliccn  oxynitride  phases  became  porous  during  the  in  situ  oxidation . 
This  type  of  reaction  can  be  seen  at  Y  in  Fig.  5  and  ir.  Fig.  6 .  The  m: r:  ho  1 1 gy  ;f  • 
phases  suggests  that  the  initial,  yttria-rich  regions  were  multiphase  or  rossit.y 
amorphous-crystalline  mixture.  The  phase  in  Fig.  6  was  identified  frerr.  its  elect.**- • 
diffraction  pattern  as  YSi02 N  or  the  K-phase  of  Fig.  1. 

The  in  situ  specimens  were  bent  and  slightly  pink  when  removed  from  th***  s'rir 
heater.  Later  TEV.  examination  shewed  that  the  oxidation  had  not  u  -nif.rmly 

ever  the  entire  thin  area  of  the  specimens.  Some  of  the  differer. ;es  re:le_t  the 
temperature  gradient  present  between  parts  of  the  specimen  tr.at  are  exp:. sed  thru-*::, 
holes  in  the  grid  heater  and  other  parts  of  the  specimen  that  were  under  the  grid 
bars.  In  the  region  shewn  in  Fig.  7  few  oxide  platelets  are  seen,  and  the  boundaries 
between  Si 3!%!*  grains  have  opened.  The  yttria  phase  visible  or.  one  edge  :f  the  grain 
was  identified  as  Y4Si2C7:.*2  or  the  J  phase.  The  cxidaticr.  rate  differences  ccul: 
reflect  also  the  inhomogenity  of  the  phase  distribution  mentioned  previously. 

Additional  in  situ  experiments  are  needed  to  verify  the  mechanism  of  the  loss  of  V.‘ 
ar.d  the  observed  reaction  differences.  Video  cr  cine"  recording  will  be  used  tc  make 
a  mere  complete  record  of  the  observations. 

The  results  of  the  in  situ  experiments  demonstrate  that:  catastrophic  low  tempera¬ 
ture  oxidation  in  yttria-silicon  material  -B-  started  at  surfaces  exposed  to  oxygen; 
voiati cation  of  a  V.*  containing  phase  created  voids;  cor.tinu :ur  films  were  produced 
over  some  grains;  yttria  rich  phases  became  porous;  am  considerable  strain  was 
ass: coated  with  the  oxidation. 
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Fig.  5.  Specimen  at  750°C  in  Torr  02 . 
The  hole  (H)  and  the  film  around  it  (Ox) 
vere  produced  in  ^30  min. 


Fig.  «.  BSi  3K4  grain  with  Si02  and  Si2K2 
platelets  produced  during  in  situ  oxida¬ 
tion.  Yttria  phase  at  Y. 


Fig.  6.  Material  -B-.  After  oxidation 
the  yttria  phase  (Y)  is  porous  and 
layered,  BSi  3^4  grains  ( B)  have  oxide 
platelets. 


Fig.  7.  Portion  of  in  situ  specimen 
vhich  was  under  the  grid,  yttria  phase 
(Y)  is  Yl4Si207N2. 
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Abstract 

TIi*»  application  of  the  coincidence  sito  lattice  (CSL)  model  to  the  diamond 
ciihc  structure  allows  a  purely  geometric  modelling  of  specific  pram  houndanes. 
Symmetric  Oil  <\SL  tilt  boundaries  with  3  were  constructed  with  oncrpies 

lower  than  the  energy  of  a  random  boundary.  Structures  without  broken  bonds 
were  found  for  rotation  angles  0<:7O  ">3\  Boundaries  within  the  ranpe  0^2G.‘>.V 
(  —  5=10)  are  considered  as  small  angle  boundaries  consisting  of  an  array  of  edge 
dislocations. 


§1.  Introduction 

There  is  considerable  interest  in  the  use  of  polycrystalline  semiconductors 
for  solar  cells  and  thin-film  devices,  and  the  investigation  of  the  electrical 
properties  of  their  grain  boundaries  is  of  increasing  importance.  One  of  the 
basic  objectives  of  this  study  is  the  determination  of  the  atomic  structure  of 
arbitrary  high-angle  grain  boundaries.  Much  progress  in  the  investigation  of 
grain  boundaries  has  already  been  achieved  in  f.c.c.  metals,  but  only  a  few 
results  were  obtained  in  covalently  bonded  crystals  (Hornstra  1959,  I960, 
Kohn  1958,  Holt  1904,  Krivanek,  Isoda  and  Kobayashi  1977).  It  is  the  pur¬ 
pose  of  this  paper  to  resume  these  investigations  in  the  diamond  cubic  structure 
from  a  more  general  approach. 

Experimental  evidence  exists  for  f.c.c.  metals  that  the  structure  of  arbitrary 
boundaries  can  be  developed  from  a  few  low-energy  special  boundaries  by 
incorporation  of  either  or  both  of  a  step-like  structure  or  misfit  dislocations. 
These  special  low-energy  boundaries  often  have  a  periodic  structure  for  specific 
misorientation  angles  of  the  two  adjacent  grains.  The  most  useful  concept 
for  the  description  of  the  geometrical  relationship  in  this  case  is  provided  by 
the  coincidence  site  lattice  (CSL)  model  (Bollmann  1970).  It  will  be  shown  in 
this  paper  that  this  model  is  applicable  to  covalently  bonded  structures  with 
the  diamond  cubic  lattice,  and  that  it  is  possible  to  develop  structures  of  specific 
(coincidence)  and  some  arbitrary  (near-coincidence)  boundaries.  Attention 
was  given  to  <01 1)  tilt  boundaries,  largely  because  a  comparison  with  models 
of  other  authors  (Hornstra  1959,  Kohn  1958)  is  possible. 
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§  2.  [01  1]  f’SL  ORIK STATIONS 

If  is  convenient  for  the  following  to  describe  the  diamond  cubic  lattice  by 
the  tet  ram  mal  unit  cell  pi  von  by  (fig.  1 )  :  a,  =  «|  1 00 1.  a,  =  n/2[Ol  1  ],  a3  =  «,'2[(>1 1 '), 
volume  I',  _  =  «3/2.  (Planes  and  directions  quoted  will  all  refer  to  the  conven¬ 
tional  diamond  f.c.c.  unit  cell — lattice  constant  a — unless  otherwise  stated.) 
When  considering  grain  boundaries  between  two  grains  with  the  diamond  cubic 


Fig.  1 


Projection  of  the  diamond  cubic  lattice  on  a  (Oil)  plane.  The  filled  and  open 
circle-  are  atom-  at  height  0  and  a  2\'2  respectively,  a,,  a2  are  basis  vectors 
of  the  direct  lattice,  b,,  b2  are  basis  vectors  of  the  CSL. 


structure,  rotation  of  one  grain  with  respect  to  the  other  around  a  common 
[Oil]  axis  can  create  a  superlattice — the  CSL — for  specific  misorientations 
which  contains  the  lattices  of  both  grains  as  sublattices.  (Oil)  is  a  symmetry 
plane  and  the  vector  rt r  therefore  becomes  coincident  with  [A,£20]trtr 

of  the  unrotated  lattice  after  a  rotation  of  0=  2  arctnn  l;2!klR  (R  =  \-'2)  (fig.  1). 
The  CSL  established  by  this  rotation  angle  9  is  tetragonal  and  can  be  described 
by  the  basis  vectors  b1  =  a[AIA-20]lctr ;  b2  =  f/[A2A,/?:!0]trtr ;  b3  =  a3  (k 2  and  k2  or 
A\2  and  k\R2  have  to  be  divided  by  common  integer  factors).  The  unit  cell  of 
the  CSL  is  similar  to  the  unit  cell  of  the  direct  lattice  with  the  ratios 

65/6,  =a,/a2=  R. 

The  volume  of  the  CSL  unit  cell  Pc  is  equal  to  2  VL  with  2  =  k22  R 2  +  A22 
(or  2  =  A  (A,2  R2  +  k22),  if  A,2  /f2-f  A-22  is  even). 

The  CSL  is  invariant  to  a  translation  of  each  sublattice  by  one  of  its  lattice 
vectors  (Pollmann  1970),  It  is  the  special  property  of  the  [Oil]  rotation  axis 
for  the  diamond  cubic  CSL  that  the  number  of  coincidence  sites  within  the  unit 
cell  of  the  CSL  can  be  increased,  in  the  eases  where  A,2  R2  +  k22  is  odd,  by  the 
translation  Ja<Ill>  of  one  lattice.  If  allowance  is  made  for  that  translation, 
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the  ratio  of  atoms  which  are  in  coincidence  is  given  by  'Z.  (100)  is  another 

symmetry  plane  of  the  lattice  perpendicular  to  the  first.  Therefore,  the  same 
('SL  can  also  be  created  by  a  rotation  of  <^=180-0. 

For  a  fixed  misorientation  of  the  two  lattices,  two  degrees  of  freedom  remain 
for  the  orientation  of  the  boundary  plane.  Periodic  boundary  structures  occur 
if  the  boundary  plane  is  a  lattice  plane  of  the  CSL.  Only  symmetric  tilt 
boundaries,  where  the  boundary  plane  contains  the  [Oil]  direction  and  is  a 
mirror  plane  of  the  CSL,  will  be  discussed  in  detail. 

♦ 

§  3.  Approximate  calculation  of  the  energy  of  a  grain  boundary 

The  wavefunctions  forming  a  covalent  bond  are  well  localized  and  it  is 
likely  that  only  a  very  limited  amount  of  relaxation  in  bond  length  or  bond 
direction  is  possible.  Therefore,  a  randomly  oriented  boundary  must  contain  a 
number  of  unsaturated  bonds.  Nevertheless,  bond  rearrangement  may  occur 
and  result  in  distorted  bonds  ;  the  approximate  energy  of  a  boundary  is 
given  by 

E  =  nEl+  X£2’(a,.),  (1) 

» 

where  Ey  is  the  mean  energy  of  a  broken  bond,  n  their  number,  and  E 2*  the 
energy  of  the  ith  distorted  bond.  Ey  is  known  for  a  variety  of  covalent 
materials  (Pauling  10(10) — for  example,  for  Si,  £’1  =  l-84eV,  and  for  Ge, 
Ey  =  I  (14  cY .  Ej  depends  on  the  distortion  angle  a,  and  can  be  approximated 
from  the  elastic  constants  of  the  crystal  under  consideration.  In  toe  diamond 
cubic  structure,  all  bonds  are  directed  along  (111,)  and  an  elastic  shear  on  a 
{111}  plane  results  in  an  angular  distortion  of  those  bonds  which  are  aligned 
perpendicular  to  this  plane.  Therefore,  at  least  for  small  distortion  angles  a(, 
the  elastic  shear  energy  per  unit  volume  V  is 

El  P  =  jp  sin2  a, 

where  p  is  the  shear  modulus  on  {1 1 1} ;  El  V  can  be  set  equal  to  the  energy  of 
the  distorted  bonds.  There  are  two  bonds  perpendicular  to  the  {111}  plane 
per  unit  cell  volume  VL  ;  therefore,  we  obtain 

^V(ai)  =  Ej  sin2  a,-.  (2) 

E2  is  equal  to  pa2/ 8  and  may  be  calculated  for  Si  or  Ge  ;  E2  =  7-2  and  6-6  cV 
respectively.  Assuming  for  the  time  being  that  eqn.  (2)  is  valid  for  all  dis¬ 
tortion  angles,  a  distorted  bond  has  to  be  considered  broken  if  E2'(oli)=E1 
which  occurs  for  a,  ft  30°  for  both  Si  and  Ge. 

The  energy  E2'( <x()  is  likely  to  increase  faster  with  angle  than  predicted 
by  eqn.  (2),  and  the  total  boundary  energy  will  therefore  be  underestimated. 
Additionally,  an  energy  term  due  to  changes  of  the  bond  length  has  been  omitted 
because  of  the  lack  of  appropriate  data.  However,  the  maximum  distortion 
of  the  bonds  across  the  boundary  (i.e.  the  bonds  connecting  the  two  grains) 
will  be  lowered  if  one  allows  the  atoms  to  relax  in  the  vicinity  of  the  boundary, 
thus  lowering  the  total  energy  again.  Therefore,  it  is  assumed  that  eqn.  (2) 
is  a  fuirlv  good  approximation  of  the  boundary  energy  at  0  K. 

Boundaries  were  formed  by  simple  geometrical  fitting  of  the  bonds  of  the  two 
lattices  across  the  boundary  plane  and,  in  general,  some  relative  translation  was 
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necessary  to  minimize  the  energv  of  tlie  distorted  bonds.  A  few  random 
boundaries  with  a  common  \U1I  axis  were  also  constructed  in  this  way,  but 
their  energy  was  very  insensitive  to  such  translation  (0=2!)  and  15(1’,  the 
table).  An  average  energy  per  unit  area  y  was  obtained  whose  values, 
~  2500  erg/ent2  for  Si  and  ~  23on  erg 'em8  for  (le,  are  lower  than  twice  the 
surface  energy  of  a  comparable!  Jlltt}  surface  of  these  materials  (Jaeeodine 
1003),  3020  erg  em8  for  Si  and  2<ioo  erg  cm-  for  Go.  Thus,  a  binding  force 
exists  between  the  two  grains.  , 

The  energy  for  t'SL  boundaries  was  in  general  lowest  if  the  boundary  was 
parallel  to  a  (100)(,si  plane,  the  boundary  with  the  shortest  periodicity  for  a 
given  — .  Occasionally,  however,  other  planes  like  pllo),,.,  or  (lln)is)  (or 
even  higher  index  planes)  had  to  be  considered.  ('SLs  up  to  £  =  33  yielded 
boundaries  with  energies  lower  than  the  average  value  of  a  random  one.  All 
calculated  values  are  listed  in  the  table  and  shall  now  be  discussed  in  detail. 


§  4.  Mooels  of  CSL  low  energy  boundaries 
The  boundaries  of  low  energy  have  common  features  and  can  be  divided 
into  two  groups. 

(a)  Rotation  angles  0°  <  9  $  70-53°.  The  few  special  boundaries  found  in 
this  range  have  £  =  l  +  2«8  («  =  l.2.  3,  4)  and  contained  only  distorted  bonds 
(figs.  2  5).  With  the  exception  of  the  £  =  3  boundary,  all  boundaries  can  be 
considered  as  an  array  of  the  same  type  of  defect—  a  five-  and  a  seven-membered 
ring  of  atoms— at  different  spacings. 


Fig.  2 

[Oil] 

V 

L— -  [  Too  ] 

[oil] 
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[011]  projection  of  the  (£  =  33)  symmetric  tilt  boundary ,  0  =  20-05°. 


f  The  random  laaindaries  had  on  average  3  1  Kinds  n*.  comparable  to  2-8  IkiihIsVi5 
for  c  {lb1}  surface. 
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Fig.  5 


[Oil]  projection  of  t!\c  (E  =  3)  symmetric  twin  boundary,  ff—IOHS1. 

The  structures  shown  here  for  the  E  =  19  and  9  boundaries  have  been  pro¬ 
posed  earlier  (Hornstra  1959).  He  constructed  a  number  of  high-angle  tilt 
boundaries  in  the  diamond  cubic  structure  as  arrays  of  edge  dislocations  and 
pointed  out  that  these  five-  and  seven-membered  ring  defects  are  edge  dis¬ 
locations,  with  a  Burgers  vector  b  =  a'2[01f|  and  a  (100)  glide  plane.  It  is 
tempting  to  suggest  that  for  all  misorientations  of  the  two  grains  0$  20-53° 
(E  =  19),  the  boundary  can  be  considered  a  small-angle  tilt  boundary  consisting 
of  these  edge  dislocations.  This  view  is  supported  by  the  observation  (Bourret 
and  Desseaux  1979)  that  this  type  of  edge  dislocation  could  be  identified  in 
small-angle  boundaries  with  0%5°  by  lattice  imaging. 

The  energy  per  unit  area  of  a  small-angle  boundary  has  been  calculated 
(Hirth  and  Lothe  1908  a)  for  the  whole  range  of  misorientations,  and  can  be 
given  by 

v-n»{5|;;c<,lh5i--l»(2!inh5l)},  (3, 

where  y0  =  ^fc/4w(l  -  v)  and  0m  =  6/27rr.  (Here  r  is  the  core  cut-off  radius  of 
the  edge  dislocation.)  This  equation  has  been  fitted  to  the  calculated  boundary 
energies  of  X  =  33  and  19  (fig.  0)  and  yielded  values  for  the  parameters  y0  = 
17  3<»»erg/cm2  and  r=  1-106.  The  first  value  does  not  agree  with  the  calculated 
y0  of  ~ 2000  erg/cm2,  but  the  cut-off  radius  is  only  somewhat  higher  than 
expected  for  an  edge  dislocation  in  the  diamond  cubic  lattice,  for  which 
7*0-36  was  proposed  (Hirth  and  Lothe  19G8b).  Equation  (3)  is  not  valid  if 
the  cores  of  the  dislocations  overlap,  which  occurs  for  9>  2(1-53°.  Consequently, 
those  boundaries  cannot  be  considered  small -angle  boundaries,  although  the  edge 
dislocation  structure  still  can  be  distinguished  in  t he  E  =  9  boundary.  Thus  it 
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Energy  (|><t  unit  area)  y  of  a  small  angle  tilt  boundary  versus  misorientation  angle 
0  (cqn.  (3)).  The  s<)iiares  are  the  calculated  C.SL  grain  boundary  energies. 
See  text  for  further  discussion. 


is  unlikely  that  t ho  energy  of  boundaries  with  6>  2(5-33°  can  be  calculated  front 
the  energy  of  their  individual  edge  dislocations.  The  high  energy  of  the  random 
boundary  (0  =  29°,  the  table)  indeed  shows  that  eqn.  (3)  is  not  valid  for 
e>  2(5-33°. 

The  first-  (-  =  3)  and  second-order  (-  =  9)  twin  boundaries  with  (100)cl51i 
plane  have  been  thoroughly  investigated  previously  (Hornstra  1939,  Kohn 
1958).  The  models  given  in  these  papers  agree  with  those  given  here;  the 
E  =  9  boundary  has  been  confirmed  by  lattice  imaging  (Krivanek  ct  al.  1977). 
It  should  be  mentioned  that  in  all  cases  studied  so  far,  relative  displacement  of 
the  two  grains  parallel  to  the  boundary  plane  was  unnecessary  to  achieve  a 
minimum  bond  distortion. 

{!>)  Rotation  angles  70-53°  <  8  ^  ISO’.  The  four  special  boundaries  found 
in  this  range  each  contained  unsaturated  bonds  (figs.  7-10).  The  E  =  3  and  9 
boundaries  are  identical  to  the  corresponding  C'SL  boundaries  with  a  rotation 
angle  of  (18O-°0)  and  a  (010)CSL  boundary  plane.  Hornstra  (1959)  and 
Kohn  (1958)  propose  slightly  different  structures  for  E  =  3,  but  their  models 
also  contain  the  same  numbfer  of  broken  bonds,  if  one  considers  that  Kohn’s 
atoms  with  five  or  three  neighbours  are  equivalent  to  one  unsaturated  bond. 
Hornstra’s  model  has  more  distorted  bonds  resulting  in  a  higher  energy  (the 
table) ;  we  therefore  believe  our  structure,  which  is  essentially  the  same  as 
Kohn’s,  is  more  likely.  In  his  investigation  of  twin  boundaries,  Kohn  pointed 
out  that  in  some  cases  one  can  achieve  a  better  geometrical  fitting  by  faceting 
the  boundary,  thus  increasing  the  number  of  coincidence  sites  in  the  boundary. 
Indeed,  the  faceted  structure  proposed  for  the  11  =  9  lateral  boundary  has  a 
lower  energy  than  the  unfaceted  boundary  plane  (see  the  table). 
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The  structure  for  2  =  11  differs  again  from  Hornstra’s  model.  He  assumes 
a  relative  translation  of  the  two  grains  along  their  common  [01 1  ]  axis  by  a'4\'2, 
which,  however,  introduces  more  distorted  bonds.  The  number  of  unsaturated 
bonds  remains  the  same  as  in  the  structure  proposed  here.  The  CSL  model 
leads  additionally  to  a  X  =  17  and  a  X  =  27  boundary  :  the  energy  of  the  former 


Fig.  8 


[Oil]  projection  of  the  (S  =  3)  lateral  twin  boundary,  9=  109'47o. 
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[Oil]  project  ion  of  the  (1=  1 1)  symmetric  tilt  boundary,  0=120-52°. 


Fig.  10 


[Oil]  projection  of  the  (2  =  0)  symmetric  tilt  boundary,  0=141-00°  (second-order 

lateral  twin  boundary). 


boundary  is  slightly  lower  than  the  energy  of  random  boundaries.  Other  CSL 
tilt  boundary  planes  than  those  presented  here  all  have  structures  and  energies 
which  are  indistinguishable  from  random  boundaries. 
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§  5.  Conclusion 

The  application  of  the  CSL  model  to  the  diamond  cubic  lattice  allows  a 
systematic  determination  of  low  energy  ‘  special  ’  boundaries.  An  array  of 
five-  and  seven-membered  rings  of  atoms  (with  or  without  a  broken  bond) 
is  the  common  feature  of  all  boundaries.  The  essential  similarity  of  the 
proposed  structures  with  those  of  other  authors  who  approached  the  problem 
differently  shows  that  the  ambiguity  involved  in  a  purely  geometric  modelling 

is  not  serious  for  this  structure.  Apart  from  the  well  established  first-  and 

1  » 

second-order  twin  boundaries,  boundaries  with  N=ll  and  IP  should  exist. 
Boundaries  within  the  range  6  ^  20-53“  (—  =  lit)  can  be  considered  as  small  angle 
boundaries,  the  five-  and  seven-membered  rings  without  broken  bonds  being 
edge  dislocations  in  the  diamond  cubic  structure. 

The  application  of  this  method  to  near-coincidence  orientations  involving 
misfit  dislocations  and  to  inclined  boundaries  is  possible  and  shall  be  discussed 
in  a  forthcoming  paper.  The  extension  of  this  approach  to  materials  crystal¬ 
lizing  in  the  structurally-related  zinc  blende  and  wurtzite  structures  is  evident. 
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ON  THE  EFFECTS  OF  DIHEDRAL  ANGLE  AND  PRESSURE 
ON  THE  DRIVING  FORCES  FOR  PORE  GROWTH  OR  SHRINKAGE 


R.  M.  Cannon 


INTRODUCTION 

At  high  temperatures,  pores  in  polycrystalline  materials  can  either 
shrink  or  grow  as  a  result  of  diffusional  transport  of  atoms  to  or  from  the 
pores.  In  principle,  two  separate  processes  can  be  considered.  In  one 
case  mass  is  exchanged  between  pores  and  grain  boundaries  (or  dislocations) 
which  leads  to  changes  in  the  apparent  density  and  bulk  dimensions  of  the 
polycrystal.  In  the  other,  mass  can  be  exchanged  between  pores  which  leads 
to  a  change  in  the  size  distribution  of  the  pores  but  not  in  the  total  pore 
volume  or  external  dimensions  of  the  polycrystal.  The  problems  of  pore 
coarsening  and  densif ication  or  cavity  growth  have  been  treated  extensively, 
although  the  problems  are  usually  treated  separately.  From  these  separate 
treatments  it  is  easily  seen  that  applied  stresses  should  affect  densifica- 
tion  more  than  coarsening,  and  that  the  occurrence  of  one  of  these  processes 
should  affect  the  driving  force  for  the  other.  However,  some  more  specific- 
coupling  between  the  various  processes  are  seen  when  they  are  considered 
together. 

It  is  generally  recognized  that  the  driving  force  for  density  changes 
depends  both  on  applied  stresses  and  on  surface  tension  forces  whereas  the 
driving  force  for  pore  coarsening  results  from  differences  in  surface  tension 
forces  at  pores  of  different  sizes.  For  isotropic  materials  attention  is 
usually  focused  on  the  pore  curvature  when  considering  these  driving  forces, 
although  it  is  recognized  that  forces  due  to  surface  tensions  can  also  be 
important.  Because  the  dominant  terms  are  usually  a  result  of  surface  curva¬ 
ture,  it  is  clear  that  for  materials  with  high  grain  boundary  energies  and 
low  dihedral  angles  the  surface  curvatures  and  driving  forces  for  densifica- 
tion  will  be  affected.  The  initial  modeling  which  was  based  only  on  curva¬ 
ture  considerations  suggested  that  low  dihedral  angles  would  reduce  the  dri¬ 
ving  force  for  densif ication,  that  pores  with  many  sides  (i.e.,  many  neigh¬ 
boring  grains)  would  have  a  lower  driving  force  for  densification,  and  that 
pores  with  many  sides  and  low  dihedral  angles  may  grow  instead  of  shrink. 
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However,  the  problem  is  incomplete  as  modeled  because  the  site  for  mass  ex¬ 
change  with  the  pere  "as  not  considered. 

In  the  present  paper  we  reconsider  the  driving  forces  for  pore  growth 
or  shrinkage  using  a  simple  two  dimensional  model.  The  results  indicate  the 
average  driving  forces  for  pore  growth  or  shrinkage  and  the  existence  of 
some  equilibrium  conditions.  The  inclusion  of  the  surface  tension  forces 
allows  a  more  specific  statement  of  the  effect  of  low  dihedral  angles  on 
densif icati on .  In  addition,  some  interesting  cases  of  coupling  between 
densif ication  and  pore  or  grain  growth  are  elaborated. 


THEORY 

1.  Basic  Considerations 

For  a  porous  polycrystal  at  low  stresses,  pore  growth  or  shrinkage  can 
result  from  defect  diffusion  between  pores  and  sites  of  lattice  discontinuity 
such  as  grain  boundaries  or  dislocations.  Within  good  lattice  the  driving 
force  for  diffusion  is  proportional  to  V(p-p  )  where  p  is  the  chemical  poten¬ 
tial  for  atoms,  and  py  is  that  for  vacancies.  If  surfaces  and  boundaries 
act  as  efficient  sources  and  sinks  for  point  defects  then  they  will  come  to 
local  equilibrium  and  provide  boundary  conditions  for  the  distribution  of 
(y-uv).  The  commonly  presumed  conditions  are  that  at  planar  grain  boundaries, 


u  -  u  =  p  =  u 


o  12 

nn 


that  across  curved  grain  boundaries. 


(1-a) 


A(u  '  Rv)  =  Au  =  <Ybfl  , 


(1-b) 


and  that  at  free  surfaces, 

p  -  uv  =  V)  =  PQ  '  °nnn  +  KysU  •  (1'c) 

where  a  is  the  normal  stress  acting  at  an  interface,  12  is  the  atomic  volume, 
<  is  the  curvature,  Yb  is  the  grain  boundary  tension,  and  ys  is  the  surface 
tension.  In  addition,  in  many  cases  there  may  be  rapid  diffusion  along  boun¬ 
daries  and  surfaces  in  response  to  Vp  as  give  by  Eqs.  (1)  assuming  point 
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defect  equilibration  is  rapid  compared  to  diffusion. 

At  equilibrium  = 0  and  u  will  be  uniform  throughout  for  a  single  com¬ 
ponent  system.  However,  equilibrium  does  not  obtain  for  many  problems  of 
interest,  and  to  treat  the  nonequilibrium,  kinetic  problems  it  is  desirable 
to  know  the  distribution  of  u-uv-  It  can  easily  be  appreciated  that  in  non- 
equilibrium  situations  the  distribution  of  p-uv  and  of  the  stresses  inside 
grains  will  depend  on  the  dominant  transport  mechanisms,  as  these  will  affect 
the  local  stresses  and  interface  curvatures.  For  a  multicomponent  system  the 
situation  becomes  more  complex  but  analogous.  These  problems  are  difficult 
to  solve  in  general,  particularly  when  considering  the  possibility  of  non¬ 
hydrostatic  stresses  or  arrays  of  pores  of  different  sizes.  Thus  it  is  use¬ 
ful  to  simplify  the  treatment  in  order  to  obtain  average  potentials  at  specif¬ 
ic  segments  of  interfaces.  The  result  of  such  a  treatment  gives  metastable 
equilibrium  conditions  where  they  exist  and  an  indication  of  the  directions 
for  and  average  driving  forces  for  densification  and  pore  growth  processes. 

When  considering  the  microstructural  evolution  of  a  porous  polycrystal 
there  are  three  processes  of  interest  to  be  considered  in  the  simplest  situa¬ 
tions.  These  are:  (1)  density  changes  which  occur  by  transport  of  atoms 
between  grain  boundaries  and  pores  (assuming  there  are  no  dislocations); 

(2)  pore  coarsening  which  involves  transport  of  mass  between  pores;  and 

(3)  grain  growth  which  results  from  the  disappearance  of  small  grains  and 
the  growth  of  large  grains  as  grain  boundaries  migrate  toward  their  centers 
of  curvature.  In  addition,  in  inhomogeneous  materials,  processes  of  grain 
shape  change  and  grain  boundary  sliding  will  also  occur.  Some  interesting 
aspects  regarding  these  processes  can  be  easily  seen  from  a  model  of  a  two 
dimensional  array  of  pores  and  grains.  Several  interesting  results  are  ob¬ 
tained  which  are  expected  to  be  general  for  three  dimensions.  In  addition, 
there  are  some  results  which  are  artifacts  of  the  simple  model,  but  which 
are  nonetheless  curiously  interesting. 

2.  Thermodynamic  Potentials 

For  a  complete  treatment  of  the  range  of  problems  of  interest,  it  is 
desirable  to  determine  the  potential  M-uv  throughout  all  regions  of  perfect 
lattice  and  the  potential  u  at  all  interfaces  and  dislocations  for  conditions 
of  non-hydrostatic  applied  stresses  and  non-uniform  internal  stresses.  The 
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chemical  potentials  can  be  found  as  derivatives  of  an  appropriately  chosen 
free  energy.  The  Gibbs  free  energy  as  usually  expressed, 

G  =  F  ♦  pV  ,  (2) 

wheTe  F  is  the  Helmholtz  free  energy,  p  the  pressure  and  V  the  system  volume, 
is  inappropriate  or  inconvenient  for  several  reasons,  but  can  be  suitably  re¬ 
defined  to  be  useful  for  the  present  problem.  If  the  stresses  are  non-hydro- 
static,  the  pV  term  does  not  adequately  describe  the  change  in  the  external 
work  potential  which  may  arise  from  additions  of  mass  or  changes  of  volume 
of  the  material.  This  has  been  discussed  by  Gibbs  and  extensively  since. 

The  reason  can  be  readily  appreciated  by  realizing  that  additions  of  mass  in 
different  ways  which  would  result  in  different  shapes  of  the  body  or  differ¬ 
ent  surface  displacements  would  cause  different  amounts  of  external  work. 
Further,  even  in  the  limiting  case  where  an  external  hydrostatic  stress  is 
applied,  the  bulk  volume  is  not  a  state  function  for  a  porous  material  and  so 
is  an  inconvenient  variable. 

For  the  present  purposes  it  is  convenient  to  define  the  system  as  the 
aggregate  of  condensed  phases,  and  it  is  sufficient  to  generalize  to  the  case 
where  the  tractions  are  not  the  same  at  all  surfaces,  but  where  along  any  par 
ticular  surface  there  is  only  a  normal  traction,  i.e.,  the  local  stresses  at 
the  surface  are 


and 


a  =  -p 
nn  r 


°nt  =  °nt '  =  0  ’ 


(3a) 

(3b) 


where  n  is  the  normal  coordinate  axis  and  t  and  t'  are  tangential  to  the  sur¬ 
face.  If  there  are  a  finite  number  of  regions  along  the  surface  where  the 
pressure  is  uniform,  i.e.  p^,  then  a  useful  free  energy  is 


G'  =  F  + 


&i4vi  ■  j TT  *  lPiavi  ■ 


(4) 


where  F  is  the  Helmholtz  free  energy  per  atom  for  the  local  composition  and 
elastic  stresses,  (o„),  and  A  is  the  atomic  volume.  The  change  in  the  work 
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potential  of  the  exterior  can  be  considered  as  a  series  of  reservoirs  at 
different  pressures.  The  total  volume  of  condensed  phases  in  the  piece  (i.e. 
the  system  volume)  is 

V  =  [av*  ,  (S) 

where  each  represents  the  volume  change  of  a  particular  reservoir  due  to 
the  system;  these  may  be  expressed  as 


AV!  =  V’ .  -  V!  , 
1  01  1 


(6) 


where  V' .  is  the  reference  volume  for  the  i*^1  reservoir  and  V!  is  the  current 
oi  th  1 

volume  for  the  i  reservoir. 

The  chemical  potential  for  any  species  can  be  defined  in  the  usual  way 
as 


Pk 


T>  P± » njj 


(7) 


The  procedure  is  useful  for  atoms  or  defects.  For  either  it  may  be  antici¬ 
pated  that  within  good  lattice  unique  values  of  the  individual  potentials  may 
not  exist.  However,  for  physically  interesting  problems  the  important  quan¬ 
tities  are  the  differences  in  chemical  potentials  between  substitutional  spe¬ 
cies,  and  these  will  be  unique.  That  is,  only  exchanges  of  species  which 
preserve  the  lattice  are  meaningful  unless  new  lattice  sites  are  created. 

The  free  energy  changes  for  these  operations  are  differences  in  potentials 
between  the  species,  and  these  are  unique  quantities.  New  lattice  sites  can 
only  be  created  at  sites  o.  lattice  discontinuity,  and  at  these  sites  the  po¬ 
tentials  of  the  individual  species  will  be  unique.  For  the  present  we  will 
assume  that  the  boundaries  and  surfaces  are  incoherent  and  good  sources  or 
sinks  for  point  defects  and  so,  at  the  interfaces,  all  species  are  at  local 
equilibrium.  (Dislocations  will  be  assumed  to  be  absent.) 

The  first  term  in  Eq.  (4)  includes  the  elastic  strain  energies  and  so 

2 

the  chemical  potentials  include  terms  of  order  a  D/2E.  For  most  sintering 
and  coarsening  problems  these  terms  are  small  and  will  be  ignored.  The  depen 


dence  of  F  on  composition  is  not  generally  negligible,  and  must  be  considered 

to  explicitly  treat  segregation  which  may  result  for  either  thermodynamic  or 

kinetic  reasons.  Presently,  we  confine  our  considerations  to  single  compon- 

f  *1 

ent  materials  or  stoichiometric  compounds!  '  Finally,  F  contains  surface 
(* ) 

'For  a  compound  with  a  narrow  range  of  stoichiometry,  under  nonequilibrium 
conditions  slight  non- stoichiometry  will  develop  primarily  near  the  inter¬ 
faces  due  to  the  differences  in  the  diffusivities.  This  will  adjust  the  po¬ 
tentials  to  cause  stoichiometric  transport.  For  our  purposes,  compounds  can 
be  considered  single  component  if  all  transfers  are  taken  to  involve  stoichio¬ 
metric  quantities. 

energy  terms  of  the  form  Jy^dA^  for  each  type  of  interface.  If  the  material 
is  sufficiently  free  of  surface  active  impurities  that  no  depletion  or  en¬ 
richment  of  surfaces  occurs,  and  if  the  interfaces  are  isotropic,  then  all 

6y.  =  0. 

J 

Thus,  the  free  energy  in  Eq.  (4)  reduces  to  that  for  a  single  component, 
isotropic,  rigid  lattice  system  of  n  atoms  (or  "molecules").  If  only  con¬ 
stant  pressure  conditions  are  considered,  then  it  is  convenient  to  use 


=  G'  -  Fp.V' .  , 

Lri  01 


so  that  G  contains  only  terms  which  are  variable,  i.e., 

G  -  nFo  -  JpjVj  .  vbAb  •  ysAs  ,  (9) 

where  A^  and  Ag  are  the  grain  boundary  and  surface  areas,  respectively. 

The  chemical  potential  for  atoms  at  a  particular  location  can  be  deter¬ 
mined  as: 


Wk  ~  <5n  T^  ”  yo 


6V. 

-fri-r- 


where  uq  =  Fq  .  Clearly,  for  a  single  component  system  the  chemical  poten¬ 
tial  will  depend  on  the  extent  to  which  work  is  done  against  a  particular  res¬ 
ervoir  and  a  particular  interfacial  area  is  changed.  In  general,  it  will  be 
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different  throughout  the  system.  If  the  surface  tension  terms  are  small,  the 
only  condition  for  equilibrium  is  if  all  the  are  equal.  If  not,  the  body 
will  have  non-hydrostatic  elastic  stresses.  These  can  be  maintained  in  mech¬ 
anical  equilibrium,  but  cause  a  chemical  potential  gradient  which  is  a  driv¬ 
ing  force  for  diffusion  which  would  lead  to  a  continuous  shape  change.  If 
the  pressure  is  not  the  same  on  all  surfaces,  the  only  way  in  which  equilib¬ 
rium  can  be  obtained,  i. e.,  y  the  same  everywhere,  is  if  the  surface  tension 
terms  become  large. 

For  the  present  paper,  the  chemical  potentials  at  different  locations 
in  polycrystalline  arrays  will  be  considered.  Attention  will  be  confined  to 
the  potentials  at  boundaries  and  surfaces.  When  these  are  equal  on  all  sur¬ 
faces,  then  no  diffusion  potentials  will  exist,  the  point  defect  concentra¬ 
tions  in  the  grain  interiors  will  become  uniform  and  at  =  0,  and  the  val¬ 
ues  of  y  in  the  interior  will  be  clearly  defined. 

3.  Geometrical  Arrays 

To  evaluate  the  driving  forces  it  is  desirable  to  consider  a  model  in 
which  both  the  size  and  the  number  of  sides  of  pores  or  grains  can  be  consid¬ 
ered  as  separate  variables.  We  consider  a  nearly  infinite  two  dimensional 
network  of  pores  and  grains.  The  network  can  be  considered  to  be  externally 
loaded  in  various  ways.  Here  we  consider  only  the  case  where  it  is  immersed 
in  a  hydrostatic  reservoir  at  pressure  pg.  To  simulate  the  situation  of 
open,  interconnected  porosity,  the  pressure  in  all  pores  would  be  equal  to  p  . 
For  the  more  interesting  case  of  closed  porosity,  the  pressures  in  different 

pores  can  be  different  from  p  and  from  each  other  and  taken  as  p  .  for  the 

.th  e  pi 

l  pore.  For  this  situation  the  free  energy  expression  in  Eq.  (10)  becomes: 

G  =  nF  -Tp.V.-pV  +  y,A^  +  y  A  .  (11) 

o  ‘•'pi  pi  re  e  b  o  s  s 

To  maintain  correspondences  with  the  three  dimensional  problem  we  can  consid¬ 
er  a  unit  depth  of  material  and  refer  to  three  dimensional  quantities  such  as 
volumes  and  pressures.  Although  the  two  dimensional  problem  is  simpler  than 
the  three  dimensional  one,  a  general  treatment  of  an  arbitrary  array  of  pores 
and  grains  is  still  formidable.  A  range  of  more  simplified  arrays  can  be 
analyzed  which  give  insight  into  the  problems  associated  with  sintering. 
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Arrays  which  are  periodic  can  be  analyzed  simply,  and  can  be  conceived  as 
part  of  an  infinitely  large  array  of  grains  and  pores.  The  simplest  arrays 
can  be  made  son.ewi  at  artificial  by  keeping  some  variables  uniform.  Tnis 
feature  leads  to  a  few  interesting  but  perhaps  unrealistic  results.  It  has 
the  advantage  that  the  effects  of  different  variables  or  free  energy  terms 
can  be  considered  separately.  With  appropriate  interpretation  most  of  the 
results  anticipated  for  realistic  arrays  can  be  seen  and  the  expected  results 
for  three  dimensional  arrays  approximated. 

Several  periodic  and  repeatable  arrays  which  fill  two  dimensional  space 
can  be  made  from  circles  of  the  same  size,  but  which  result  in  pores  with 
different  sizes  or  numbers  of  sides.  Kith  reference  to  the  sintering  problem 
it  is  easy  to  conceive  of  these  arrays  as  resulting  from  the  packing  of  cir¬ 
cles  (or  rods)  in  the  specified  arrangements  and  sintering  beyond  the  initial 
stage.  Qualitatively,  there  are  several  different  types  of  arrays  which  dif¬ 
fer  in  the  degree  to  which  the  local  geometry  can  be  easily  characterized  or 
specified.  There  are  three  arrays  in  which  the  grains  are  identical  and  all 
of  the  pores  are  identical  and  for  which  the  grains  and  pores  are  symmetrical. 
These  can  be  obtained  from  hexagonal  packing  of  circles,  square  packing  of 
circles,  and  hexagonal  packing  of  circles  with  1/3  of  the  circles  missing. 
Figure  1.  These  are  characterized  as  3- (6+3) ,  4- (4  +  4),  and  6- (3+3) ,  respec¬ 
tively,  where  the  first  number  is  the  number  of  sides  per  pore  and  the  number 
of  sides  per  grain  is  given  in  parentheses  as  the  sum  of  the  number  of  grain 
boundaries  and  pore  faces.  The  second  class  of  arrays  considered  here  are 
those  in  which  n'  circles  are  packed  in  repeatable  cells  of  n'+l  sites.  Fig.  2. 
These  give  a  combination  of  one  six-sided  pore  and  several  three-sided  pores 
per  cell.  These  arrays  are  also  amenable  to  considerable  analysis  because 
both  types  of  pores  are  symmetrical  and  the  geometry  can  be  plausibly  speci¬ 
fied  and  characterized.  Arrays  in  which  the  pores  are  asymmetric  or  have 
more  than  one  type  of  facet  are  of  interest,  but  can  only  be  treated  more 
approximately.  These  include  arrays  of  5- sided  pores  with  equal  numbers  of 
7-sided  grains  and  packing  accidents  involving  clusters  of  missing  particles, 
i.e.  divacancies,  trivacancies,  etc.,  up  to  large  holes. 

4.  Triple  Point  Equilibria 

For  two  dimensional  networks,  only  junctions  between  three  interfaces 


are  stable.  Although  junctions  of  four  or  perhaps  more  interfaces  would  form 
occasionally  under  non-equilibrium  conditions,  they  are  sufficiently  unstable 
that  they  would  decompose  rapidly.  Thus,  we  only  consider  networks  with 
three  interface  intersections.  The  exception  to  this  which  is  of  some  inter¬ 
est  occurs  when  large  pores  disappear  which  would  lead  to  an  instability. 

At  three  interface  junctions,  local  force  balance  or  energy  minimization 
dictates  that  the  interfaces  meet  at  the  equilibrium  dihedral  angles.  For 
three  grain  boundary  junctions,  this  requires  120°  angles  between  grain  boun¬ 
daries.  For  grain  boundary-pore  intersections,  the  equilibrium  angle  depends 
on  the  ratio  of  'V,j)/ys  as 


cos(4>/2)  =  Tb/2ys  .  (12) 

(The  angle  4>  is  defined  in  Fig.  3.) 

5.  Pore  Geometry 

Under  non- equilibrium  conditions,  the  chemical  potentials  will  vary  con¬ 
tinuously  along  all  interfaces  (assuming  they  are  good  point  defect  sources), 
and  the  chemical  potentials  must  be  continuous  at  triple  points.  Thus,  from 
Eqs.  (lc)  and  (la)  it  can  be  seen  that  the  radius  of  curvature  of  the  pores 
will  vary  continuously  as  will  the  normal  stresses  acting  on  boundaries.  The 
actual  pore  shapes  and  stress  distributions  along  boundaries  will  depend  on 
the  dominant  transport  mechanism.  However,  at  triple  points  very  little  atom¬ 
ic  transport  is  needed  to  establish  the  equilibrium  dihedral  angle.  Thus,  it 
is  both  energetically  favorable  and  kinetically  possible  that  the  equilibrium 
dihedral  angles  will  establish  quickly  even  at  the  expense  of  inducing  higher 
curvature  in  adjoining  segments  of  some  interfaces. 

At  equilibrium,  uniformity  of  the  chemical  potential  along  the  pore  sur¬ 
faces  requires  that  each  pore  segment  be  a  circular  arc.  Similarly,  the 
stresses  along  the  boundary  must  be  uniform.  These  conditions  follow  direct¬ 
ly  from  calculations  of  the  requirements  for  complete  equilibrium.  To  obtain 
annroximate  driving  forces  for  transport  in  non- equilibrium  conditions  req¬ 
uires  rigorous  kinetic  calculations  or  some  simplifications.  To  do  the  lat¬ 
ter  here  we  assume  that  the  pores  are  bounded  by  segments  of  circles  which 
meet  at  the  equilibrium  dihedral  angle  at  the  grain  junctions.  With  those 
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assumptions,  average  chemical  potentials  can  be  calculated  for  pore  surfaces 
or  grain  boundary  segments  from  Eq.  (11).  Although  tl 3  normal  grain  boundary 
stresses  need  not  be  calculated  explicitly,  only  average  potentials  are  ob¬ 
tained  which  are  equivalent  to  those  which  would  result  assuming  the  boundary 
segment  were  uniformly  stressed. 

For  symmetrical,  equilibrium  pores,  the  sign  and  magnitude  of  the  radius 
of  curvature  can  be  explicitly  determined  and  will  depend  on  the  pore  area, 
the  dihedral  angle,  and  the  number  of  sides  of  the  pore,  N.  The  pertinent 
geometrical  terms  are  defined  in  Figure  3.  The  chemical  potentials  are  most 
conveniently  described  in  terms  of  the  radius  of  curvature,  R.  For  large  di¬ 
hedral  angles  the  pore  surfaces  will  be  concave  and  we  take  R  <  0  for  this 
situation.  For  small  dihedral  angles  the  pore  surfaces  will  be  convex,  for 
which  R  >  0.  Obviously,  for  some  dihedral  angle,  the  pores  will  have 

straight  sides,  i.e.  R  ■+•  ”.  This  critical  angle  will  be  larger  for  pores 
with  many  sides,  i.e., 


*c  =  rrCl  -  2/N)  . 


(13) 


For  a  given  N,  the  dependence  of  R  on  pore  size  and  4>  can  be  seen  from: 


R  = 


-S 

2  sin  Aij>/2  ’ 


(14) 


where  A<fi  =  The  size  of  the  pore  is  more  conveniently  visualized  in 

terms  of  r,  the  radius  of  the  circle  which  circumscribes  the  pore.  However, 

the  relationships  between  R  and  the  pore  area  A,  or  r,  are  more  cumbersome. 

We  note  that  as  N  increases  or  4>  decreases  the  ratio  r/R  increases  algebraic- 

2 

ally  (note  that  -1  <  r/R)  and  that  the  ratio  A/  nr  decreases  as  N  or  ♦  de¬ 
crease. 

As  will  become  apparent,  in  some  cases  it  is  useful  to  describe  the 
results  either  in  terms  of  the  fractional  porosity  or  in  terms  of  the  initial 
configuration  of  circles  used  to  construct  the  arrays.  For  the  arrays  con¬ 
sidered  which  are  made  up  initially  of  n*  or  (n'+l)  circles  of  diameter  Dq  on 
(n'+l)  sites  the  porosity,  v^,  can  be  expressed  as 

n  .  /r.\2 

vf  =  £  (n'  +  l)  (T/  fi^  ’ 


(15a) 
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or,  unless  R  =  ®,  as 


n  .  /R.\2 

vf  =  ^  (n’+l)\Tj  ®i  ^  ’ 


(15b) 


where  there  are  nci  pores  of  a  particular  type  (i.e.  N)  and  size  in  the  cell; 
the  geometrical  factors  depend  on  N  and  the  packing  configuration.  The  grain 
size,  D,  is  conveniently  defined  in  terms  of  the  diameter  of  the  sphere  that 
could  occupy  the  site,  Figure  4.  This  is  related  to  Dq  by 


D  = 


(16) 


where  v^q  is  the  porosity  of  the  initial  array  which  is  the  maximum  possible 
for  that  particular  packing  configuration.  Also  note  that  A^,  the  initial 
pore  area,  depends  strongly  on  N  and  is  the  maximum  possible  for  a  pore  of  N 
sides. 


RESULTS 

Average  driving  forces  for  different  processes  can  be  determined  by  con¬ 
sidering  the  variation  in  G  as  given  in  Eq.  (11)  for  various  situations  or 
arrays.  From  Eq.  (11)  at  constant  p^  and  pg  , 

6G  =  F6n  -  Ipp.6V.  -  p^  ♦  Y^  ♦  Y^  .  (17) 

By  determining  6G  subject  to  appropriate  constraints,  chemical  potentials  or 
chemical  potential  differences  between  two  locations  can  be  determined.  It 
is  useful  to  first  consider  mass  additions  to  the  polycrystalline  arrays  in 
order  to  indicate  the  origin  of  specific  driving  forces.  Then  mass  exchanges 
between  locations  will  be  considered  which  represent  densif ication,  pore 
growth,  or  coarsening  processes  of  interest.  In  several  cases  equilibrium 
conditions  are  obtained  for  some  processes,  but  for  most  cases  these  are  only 
metastable  or  quasi-equilibrium  conditions. 

1.  Simple  Mass  Additions 

First  consider  an  addition  of  mass  from  an  external  source  to  a  particu¬ 
lar  pore  surface.  For  this  process 
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6n  =  -SVp./fi  ,  (18) 

and  for  this  process  all  other  pores  and  the  external  reservoir  are  left  un¬ 
affected,  i .  e. , 


6V  .  =  SV  =0 
pj  e 


(19) 


For  all  of  the  configurations  considered,  it  can  be  shown  that,  if  the  di¬ 
hedral  angle  is  the  equilibrium  value  given  by  Eq.  (12),  then 


<$n 


*  Y, 


6A 


6n 


V 

R. 

1 


(20) 


and  from  Eqs.  (17-20), 


p  -  p  =  -= —  +  p  .n  . 

P  o  Ri  *pi 


(21) 


As  can  be  seen  in  Figure  3,  regardless  of  the  sign  or  value  of  R,  it  is  true 
that  fiA^/finiO,  6As/6n<0,  and  <S  ( R [ /6n  <  0.  Thus,  at  p  ^  =  0  the  magnitude  of 
Pp-PQ  will  increase  as  4>  is  different  than  <j>c>  and  will  be  greater  for  small¬ 
er  pores;  however,  the  sign  of  pp-pQ  will  depend  on  whether  <f>  is  large  or 
small  relative  to  4>c>  It  is  worth  emphasizing  that  addition  of  mass  to  a 
particular  pore  does  not  result  in  changes  in  the  conditions  elsewhere.  Thus 
even  if  large  additions,  An,  are  made  to  the  pore,  p^  depends  only  on  R  and 
Pp  for  the  particular  pore. 

Alternatively,  consider  an  addition  of  mass  to  the  grain  boundary  re¬ 
gion.  In  the  rigid  lattice  approximation  even  infinitesimal  amounts  of  mass 
cannot  be  added  to  a  local  region  of  grain  boundary  or  even  around  an  entire 
grain  because  prohibitively  high  stresses  would  be  developed.  There  are  two 
alternative  ways  to  proceed.  One  is  to  consider  the  additions  spread  along 
the  boundaries  in  such  a  way  that  no  incompatibilities  develop.  The  simplest 
way  is  to  consider  a  uniform  addition  to  all  boundaries  as  shown  in  Figure  4. 
The  effect  of  the  addition  on  the  pore  surface  or  volume  is  somewhat  arbi¬ 
trary.  If  it  is  prescribed  that  the  addition  occurs  such  that 


6V  .  -  0 
Pi 


(22) 
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then 

6As.  =  0  ,  (23) 

and 

6n  =  -fiV£/n  .  (24) 


(25) 


(26) 


The  second  term  is  expected.  The  first  terra  has  generally  not  been  consider¬ 
ed,  but  the  schematic  in  Figure  4b  clearly  indicates  that  if  mass  is  added 

to  the  system  such  that  6V  <  0  and  the  number  of  grains  in  the  array  is  con- 

c 

stant  then  6 A,  >  0. 

b 

Alternatively,  a  sufficiently  small  layer  of  mass  may  be  added  locally 
to  a  grain  boundary  region  if  it  can  be  accommodated  elastically.  The  change 
in  external  volume  will  be  that  given  in  Eq.  (24)  which  would  lead  to  the 
second  term  in  Eq.  (26)  for  The  correspondence  with  the  first  term 

in  Eq.  (26)  is  harder  to  establish.  If  the  added  mass  is  considered  as  an 
interstitial  dislocation  loop  along  a  region  of  boundary,  then  there  are  two 
terms  which  could  be  identified  with  an  increase  in  A^.  Depending  upon  where 
in  the  net  the  loop  terminates,  the  edge  may  be  associated  with  an  increase 
in  which  formally  leads  to  a  term  involving  y^SA^.  In  addition,  other 
boundaries  in  the  network  will  be  stretched  elastically  as  part  of  the  elas¬ 
tic  accommodation  of  the  loop.  However,  the  work  done  to  stretch  the  inter¬ 
face  elastically  is  not  generally  the  same  as  y^,  which  results  from  increas¬ 
ing  A^  by  adding  more  atomic  sites  and  leaving  the  surface  structure  statisti 
cally  unchanged.  Further  analysis  to  establish  the  conditions  under  which 
this  process  gives  the  potential  in  Eq.  (26)  is  beyond  the  present  scope,  but 
we  make  several  observations.  The  local  additions  of  mass  will  lead  to  sig¬ 
nificant  stress  redistribution  if  in  becomes  appreciable  which,  of  course, 
must  be  properly  treated  in  transient  processes.  In  many  cases  this  redistri 
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bution  is  slow  enough  that  Eq.  (la)  will  still  describe  in  terms  of  the 

local  strerses  However,  if  the  minimum  admissible  6n  is  sufficiently  large 

that  the  stresses  would  be  appreciably  affected,  then  neither  p  £1  nor  o  fi 

e  nn 

may  adequately  describe  the  external  work  and  significant  changes  in  the 
2 

Jo  dV/2E.  This  would  seem  most  likely  to  be  a  problem  in  describing  problems 
of  relaxation  of  high  local  stress  concentrations  or  pore  nucleation.  Most 
of  the  sintering  and  coarsening  problems  of  interest  are  quasi-steady  prob¬ 
lems  in  which  the  expression  for  in  Eq.  (26)  gives  a  correct  view  if 

the  volume  over  which  mass  is  spread  is  properly  considered.  In  these  cases 
the  additional  mass  clearly  leads  to  an  increase  in  the  number  of  grain  boun¬ 
dary  sites  of  excess  free  energy 


2.  Densification 

Densification  results  from  transfer  .of  mass  from  grain  boundaries  to 
pore  surfaces  and  results  in  a  change  in  the  external  dimensions.  The  poten¬ 
tial  difference  between  the  pores  and  grain  boundaries  can  be  determined  from 


U 


P 


6G  __  6G 

6n  6n. 

P  1 


(27) 


This  is  subject  to  the  constraint  that  the  system  volume  does  not  change, 
6V  =  0,  which  leads  to 


[6V  .  =  -<5V  .  (2< 

‘-pi  e 

If  exchange  between  all  the  boundaries  in  a  cell  and  all  the  pores  of  one 
particular  type  in  the  cell  is  considered,  then 


<5n  .  =  -6V  /a  =  -5n.  =  SV./fi  , 
pi  pi  b  b 


and  from  Eqs.  (17),  (27)  and  (29), 


V 


Pi 


(29) 


(30) 


A  set  of  relations  like  Eq.  (30)  exists  for  each  type  of  pore 
Jb 


Obviously, 


Upj-wK  can  be  positive,  negative  or  zero,  which  describes  the  tendency  of 
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the  pore  to  grow,  shrink,  or  be  at  equilibrium  with  respect  to  the  grain  boun¬ 
dary. 

It  is  possible  for  a  particular  pore  or  set  of  pores  to  be  at  equilibrium 
with  respect  to  the  grain  boundaries,  i.e.  p  =  p^,  when  the  bracketed  term 
in  Eq.  (30)  vanishes.  Several  specific  situations  are  worth  noting. 

If  <p  -  n,  then  R  =  -r  and 


-  w.  = 
P  b 


(yb  1 s 


pP  r  • 


Obviously  p^-p^  <  0  unless  there  is  an  internal  pressure,  p^  >  0,  or  an  exter¬ 
nal  tension,  i.e.  -p^  =  o  >  0  .  Unless  p^  >  0  or  a  is  sufficiently  large, 
sintering  to  full  density  is  the  equilibrium  state. 

If  <f>  <  <J>  ,  then  R  >  0  and  p  =  p,  =  0  is  possible  even  with  p  =  p  =0. 
c  p  b  1  re  rp 

From  Eq.  (30)  the  equilibrium  condition  for  the  stress-free  case  is  when 


and  more  generally  when 

D  2cos ($/ 2)  -  (pp-Pe)D/Ys  ’  1 

so  if  the  pores  are  large,  R  >  R*  and  they  tend  to  shrink.  If  the  pores  are 
small,  R  <  R*  and  they  tend  to  grow.  At  R  =  R*  the  pore  is  at  a  metastable 
equilibrium.  The  tendency  to  growth  or  shrinkage  depends  both  on  the  pore 
curvature  which  depends  on  and  N  and  on  the  grain  size.  This  is  shown  in 
Figure  5  for  the  stress-free  case.  Interestingly,  for  <#>  >  0  some  shrinkage 
is  indicated  for  all  pores. 

The  question  of  what  occurs  when  N  is  very  large.  Figure  6,  is  of  some 
interest  as  R  >  0  even  for  quite  high  By  considering  very  large  pores  it 
can  be  shown  that 


-  u. 
P  b 


where  r  is  the  radius  of  the  pore.  Thus  even  large  pores  tend  to  shrink 
initially. 
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From  Eq.  (15)  it  is  clear  that  a  particular  value  of  R*/D  corresponds 
to  a  particular  porosity.  Thus  Eqs.  (32)  and  (33)  indicate  tnat  whe  i  $  <  <J> 
there  is  an  "equilibrium"  porosity  which  depends  on  the  number  of  sides  for 
the  pore  and  on  the  internal  pressure  and  applied  stress.  The  effect  of  <f> 
and  N  on  the  equilibrium  porosity,  v£,  at  =  pg =  0  is  shown  in  Figure  7. 
The  effect  of  pressure  can  be  seen  from 


v*(o  +  p)  1 

v*(0)  [1  -  (o  +  pp)D/Yb)2 


(35) 


which  is  plotted  in  Figure  8. 

An  interesting  result  on  bloating  in  the  presence  of  an  inert,  entrapped 
gas  can  be  anticipated  from  Eq.  (34)  or  Figure  8,  which  show  that  when  the 
equilibrium  porosity  results  from  p^  >  0,  the  relation  depends  on  grain  size. 
The  effect  can  be  seen  clearly  by  assuming  the  gas  is  ideal,  in  which  case 


n  kT  n  kT 


(36) 


As  densification  occurs  then  p  will  increase  until  p 

P 

tion  will  stop.  This  will  occur  when 


P  Mb 


0  and  densifica- 


P 


P 


n  kT 
=  -S -- 
vf  y 


(37) 


Further  annealing  will  result  in  grain  growth  and  pore  growth.  If  both  R 
and  D  increase,  the  equality  in  Eq.  (36)  cannot  be  maintained  if  the  porosity 
remains  constant,  i.e.,  if  R/D  remains  constant.  Thus,  as  R  and  D  increase, 
Vj.  must  increase  which  requires  that  R  increases  faster  than  D.  This  will 
lead  to  sintering  behavior  in  which  the  density  goes  through  a  maximum,  as 
shown  schematically  in  Figure  9. 


SUMMARY 

The  thermodynamics  for  pore  growth  or  shrinkage  have  been  reformulated. 
The  chemical  potentials  which  describe  pore  growth  or  shrinkage  have  terms 
which  depend  on  the  pore  curvature,  the  line  tension  of  the  grain  boundaries, 
the  applied  stress  and  internal  pressure  in  the  pores.  The  effect  of  pore 
curvature  may  be  positive  or  negative  depending  on  the  dihedral  angle.  The 
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additional  grain  boundary  tension  term  always  favors  densif ication. 

The  results  using  a  simple  model  indicate  that  at  low  dihedral  angle 
(<j>  <  60°  in  two  dimensions,  $  <  70.5°  in  three  dimensions)  there  is  an  equi¬ 
librium  porosity  in  a  polycrystalline  material;  for  such  a  material  there  is 
no  barrier  to  nucleation  of  cavities  in  dense  material. 

At  intermediate  dihedral  angles,  pores  with  few  sides  would  shrink  and 
disappear,  whereas  those  with  many  sides  would  initially  shrink  but  would 
reach  an  equilibrium  size  at  which  densification  would  stop.  In  this  case 
grain  growth  without  growth  of  the  large  pores  would  reduce  the  number  of 
sides  for  the  pore  and  promote  further  densification.  For  such  materials 
there  would  be  a  small  barrier  to  nucleation. 


Figure  1.  Hexagonal  array  of  grains  with  only  2/3  of  possible  grain  sites 
I  occupied.  This  can  be  regarded  as  a  limit  of  the  series  of  arrays  with  n'- 

sided  pores  and  (12-n' )-sided  grains.  Alternatively,  it  is  the  limit  of 
stable  cells  with  n  grains  occupying  n+1  sites.  The  evolution  of  pore  and 
grain  shapes  with  sintering  can  be  seen  in  the  center  pore. 

I 
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Figure  2.  Several  cells  for  arrays  of  n'  grains  occupying  n'+l  sites. 


(a) 


(b) 


Figure  4.  Schematic  illustration  of  the  effect  of  adding  mass  to  the  grain 
boundaries  where  it  is  distributed  uinformly  on  all  grain  boundaries  within 
the  cell.  When  all  6V  =  0,  the  method  of  addition  in  (b)  is  formally  equiv¬ 
alent  to  mass  addition^at  the  boundaries  and  better  illustrates  the  resulting 
changes,  although  it  is  physically  impossible. 
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Figure  8.  Ratio  of  the  equilibrium  porosity  to  the  maximum  porosity  as  a  function  of  applied  and 
internal  pressures  for  various  dihedral  angles. 
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